Scaling of excitations in dimerized and frustrated spin-1/2 chains 
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We study the finite-size behavior of the low- lying excitations of spin-1/2 Heisenberg chains with 
dimerization and next-to-nearest neighbors interaction, J2. The numerical analysis, performed 
using density-matrix renormalization group, confirms previous exact diagonalization results and 
shows that, for different values of the dimerization parameter 8, the elementary triplet and singlet 
excitations present a clear scaling behavior in a wide range of £ — L/£ (where L is the length of the 
chain and £ is the correlation length). At J2 = Ji c , where no logarithmic corrections are present, 
we compare the numerical results with finite-size predictions for the sine-Gordon model obtained 
using Luscher's theory. For small 8 we find a very good agreement for £ > 4 or 7 depending on the 
excitation considered. 
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Spin-1/2 chains with frustration and explicit dimer- 
ization have attracted a lot of attention for their experi- 
mental relevance to spin-Peierls materials! as well as for 
their interesting theoretical properties, related, for in- 
stance, to the presence of two independent mechanisms 
for spin-gap generation. The difficulties in studying these 
systems are a consequence of the lack of quantitative an- 
alytical methods to describe their full phase diagram and 
of different types of limitations in the currently available 
numerical methods. On the one hand, with exact di- 
agonalization one can treat at most few dozens of sites, 
which is inadequate when the correlation length becomes 
large. On the other, using the density-matrix renormal- 
ization group (DMRG) technique^ it is possible to reach 
hundreds of sites, but the control on the quantum num- 
bers of the excitations is much more complicated. In 
both cases it is difficult to reconstruct the infinite-size 
spectrum in absence of a theory that describes the finite- 
size effects, especially because some excitations (like the 
singlet - sec below) arc known to have a non-monotonic 
dependence on the chain length L (Ref. 0). In this pa- 
per we present a detailed finite-size analysis of DMRG 
results for the low-energy excitations of a spin-1/2 chain 
with frustration and explicit dimerization, with the aim 
of getting a deeper theoretical understanding of the cor- 
rections due to finite L. In particular, along a specific line 
in parameter space, where the spin chain is described by a 
massive intcgrablc quantum field theory (QFT), we com- 
pare the numerical results with analytical predictions for 
the finite-size scaling based on Luscher's theory^. 

The Hamiltonian of a frustrated spin-1/2 chain with 
explicit dimerization has the standard form 

L L 

= J ^ [(1 + (-)'*)] Si ■ S - t +i + S » ' S *+2 M 

i=l t=l 

(periodic boundary conditions are assumed in the fol- 
lowing). For <5 = one identifies two regimes. For 
J2 < 'he — 0.24J the model is gapless and described, 
at low energies, by the SU(2) Wess-Zumino-Witten- 



Novikov (WZWN) conformal field theory at level k = 
1 (SU (2)i), with a marginally irrelevant perturbation. 
When J2 > Ji c the marginal perturbation changes sign 
and induces a gap in the spectrum. In terms of operators 
of the WZWN theory the dimerization term is propor- 
tional to the trace of the SU(2) matrix field g. Then, 
in the low energy limit, Q is described by the following 
Hamiltonian density^ 



H = H 



SU{2) 



- 7 J • J + 77 Trg, 



(2) 



where Ti 



SU(2)t 



If- V. J 



J • J : , J and J 



are 

the chiral SU(2) currents and satisfy the level k — 1 
Kac-Moody algebra. The coupling constants are linearly 
related to J2 and 5: 7 cx (1 — J 2 1 J 2c) and 77 oc S. It 
should always be kept in mind that, roughly speaking, 
the field theory description is valid if the masses are much 
smaller then the bandwidth, that implies 5 <C 1 or. for 
6 = 0, (J 2 - J 2c )/J<l. 

The theory @ is integrable when one of the two cou- 
pling constants vanish. For 77 = excitations are massive 
(massless) spinons for 7 < (7 > ). When 77 ^ 0, the 
strongly relevant perturbation Trg introduces a confining 
interaction between the spinons that are no longer funda- 
mental excitations of For 7 = the model is again 
integrable and equivalent to the sine-Gordon (SG) model, 
C = l/2(d fi Lp) 2 - .9i/(2tt 2 w) cos(/fy), at 1 = 2tt, and the 
spectrum contains triplet and sing let excitations^. The 
ratio between their masses, R = m s /mt, is known exactly 
to be R = 

The behavior of R for finite values of 7 (and 77 7^ 0) is 
difficult to address analytically even in the regime of va- 
lidity of the field theory description, because (J2J, like its 
lattice counterpart |JIJ, is not integrable. This problem 
was studied using exact diagonalization in Ref. y, where 
it was found that the singlet is stable for any J2 < J^c 
and R interpolates between V3 at J2 = J2c and R = 2 
at J 2 = 0. These results are reproduced by our DMRG 
analysis, although we do not report them here (see how- 
ever Ref. for the extension to J2 < 0) . Some analyt- 
ical predictions can be made in the vicinity of the two 
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integrable points, where a pcrturbative analysis can be 
carried or&S. 

For the analysis of the above problem as well as many 
others it is essential to have a theory for the finite-size 
scaling that allows one to extrapolate numerical results. 
For massive integrable QFTs it is possible to express the 
leading finite-size corrections to the spectrum in terms 
of their exact scattering data^£. As we recalled, there 
are two points for which (J2J) is integrable. In this pa- 
per we consider the integrable point 7 = and r\ 7^ 0, 
rather than 77 = and 7 < 0, because the former case 
presents bound states and is physically more interesting. 
Numerically, this is also the easiest scenario to investigate 
because in the other case the gap opens very weakly^. 
Hence, in the following we will compare DMRG data for 
the lattice Hamiltonian with J2 = J2 C and 6 -C 1 with 
predictions coming from the SG theory at /3 2 = 2tt. 

Let us start presenting the numerical results. Before 
studying the finite size-effects we extract from the nu- 
merical data some important non-universal pre-factors 
- a problem recently addressed in Ref. It is well 
known that, in absence of logarithmic corrections, the 
triplet mass scales as m t — A t S 2 ^ 3 , with A t being a 
non-universal amplitude. The energy gaps, computed 
with DMRG, are linearly related to the masses through 
a velocity pre- factor that we express as v = C v Jn/2: 
AE t = vm t (the lattice spacing is formally set to one). 
For J2 = it is known that C v = 1, but for Ji = Ji c the 
constant may assume a different value. Another unknown 
constant connects S with the coefficient of the relevant 
term in the SG model: g x = 6J5(n/2) 1 / 4 Cs (Cs = 1 for 
J2 = 0i£). Along the same lines of Ref. Il3l we can extract 
these constants from the pre-factors of the triplet gap and 
ground-state energy density, that we write as AE t / J = 
1.723K t S 2 / 3 and [e GS (S)-e GS (0)]/J = -0.2728X GS <5 4/3 , 
respectively. The relation between the constants is K t = 
C v {C s /C v ) 2/3 and K GS = C v (Cs/C^ 3 . We have com- 
puted these quantities for some values of 5 using the 
finite-system DMRG algorithm^ with three "zips" in a 
range of L from 10 to 100. The results are presented in 
TablcHJ and imply C v ~ 0.81 and Cs — 1.53 at J2 = J20 
to be compared with Cs = C v = 1 at J 2 = 0. 

Let us now come to the central part of the paper 
by examining the finite-size features of the low-energy 
spectrum. In the framework of the DMRG this is done 
by building the superblock density matrix as a mixture 
(with equal weights) of the matrices associated with a 
prescribed number of target states. By means of the 
thick-restart variant of the Lanczos method 1 ^ wc are 
able to target up to ten excited states in each sector of 
^tot = Si^j which is the only good quantum num- 
ber that we have exploited. First of all from the so- 
called 1-particle (IP) excitations of the triplet we can 
get an independent estimate of the velocity. In fact, 
these excitations for small momenta q n = 2-Kn/ L (with 
n = ±1, ±2, . . .) have energy E n = E GS + y/ AE 2 + v 2 ql . 
So we have checked that L^J AE% — AE 2 /2tt actually 
scale as nv\ and we have reported the values of v\ in 



6 \AEt/J\Kt \K G s I vi/J \Lmi n (S) 
0.01 1 0.09875 1 1.235 1 1.966 1 1.1921 50 
0.02 0.15705 1.237 1.948 1.1822 ± 0.0010 30 
0.10 0.46233 1.245 1.800 1.1425 ± 0.0003 10 



TABLE I: Triplet gaps, constants K t , Kas and estimates 
of the "1-particle velocity" vi for different values of S and 
J2 = J2 C (see text after Eq. ffllUB for the determination of 
AE t / J.). Typically M = 256 optimized DMRG states have 
been retained. When AEt < 0.1 we have improved the preci- 
sion using M — 405. The last column contains the (approxi- 
mate) values at which the singlet gaps attain a minimum. 

Table [I] Unfortunately, these states with nonzero mo- 
mentum represent also the main limitation that we have 
encountered in the computation of the singlet gap AE S 
at large values of 8. This is because all the excited states 
that we can target in S* ot = 0, say No, are eventually 
exhausted by the triplet and its IP excitations. In prac- 
tice, the singlet "disappears" from the set of DMRG lev- 
els for L > (N - 2)irv 1 /AE t ^R 2 - 1, where R is the 
mass ratioi^. In addition, unlike the triplet, the singlet 
level has a non-monotonic dependence on L. Typically 
AE S (L) decreases with increasing L up to a characteris- 
tic value L m i n (S), beyond which it starts to increase and 
converges to the asymptotic limit from below. An exam- 
ple of such a feature is plotted in Fig. ^ for 6 = 0.02 
and J2 = J2c- As listed in Table [l]t;h e values of L m i n 
increase for 5 — > so that there exists a minimum value 
of 6 below which L min becomes larger than the maximum 
reliable size that we can handle with the DMRG. 

There are other cases in the literature for which numer- 
ical results show that some states in the spectrum exhibit 
a non-monotonic dependence on L (see for instance Refs. 
[T^l . In such situations it is particularly difficult to get a 
reliable estimate of the infinite-size value. In some cases 
it is possible to push L beyond the minimum and even 
beyond the flex where the finite-size data start to satu- 
rate. Whenever this is not the case the construction of a 
suitable fitting function becomes essential. In Rcf.slsHT^. 
some empirical scaling functions were suggested that pro- 
vide a good fit of the numerical data. Nevertheless these 
expressions were not obtained from a consistent theo- 
retical construction and contain various fitting param- 
eters. In the following we will construct the finite-size 
corrections using the QFT description. In particular we 
will compute rigorously the dominant contribution to the 
scaling functions for large L. These expressions contain 
no fitting parameters. 

As already recalled for J2 — J2 C and S <g; 1 the spin 
chain is described by the integrable SG model at 1 = 
2ir. The factorized scattering theory explicitly shows the 
hidden SU(2) symmetry of the model at this pointi. The 
spectrum consists of a soliton s and an anti-soliton s 
which are degenerate with a breather state b\. Moreover, 
all these particles have the same S-matrbAi£ 

q ia\ q ta\ sinhfl + i sin tt/3 _ 

i>a,b{p) = do{V) = — — tx> a,b = s,s,b 1 . (3) 

smhO — ism 71-/3 
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FIG. 1: Finite-size gaps vs 1/L for 8 = 0.02 and J2 = Jic- 
triplet (open squares), singlet (full circles), lowest IP triplet 
excitations with momenta different from and 7r (coinciding 
up and down triangles). Following the singlet at large L re- 
quires to target various excited states within S* oi = (No — 6 
in this example, where M = 256), due to the crossing with IP 
states. The singlet curve has a minimum at L m in(0.02) ~ 30. 

The rapidity 9 parameterizes the relativistic dispersion 
relations: e a = m a cosh(9, p a = m a smh9. In addition, 
there is another breather state &2 of higher mass. The 
additional S-matriccs aroiiS 

Sa M (d) = S Q (6 + m/6)S (6 - in/6) , a = s, s, 61 ; (4) 
S b2 , b M = [S o (0)} 3 . (5) 

The excitations can be then organized into a triplet 
(s, s, bi) of mass m t and a singlet of mass m s = \f3mt. 

For what follows, it is worth recalling that bound 
states, like the breathers, are associated to poles in the 
S-matrix. In particular if 9 = iu^ c is a pole in the scat- 
tering process of the particle b and c the mass of the 
bound state is given by: m 2 , = m\ + to 2 . + 2m b m c cos u% c . 
For instance, the simple poles in the S-matrix S StS (9) at 
9 = 2iir/ 3 and 8 = iir/3 correspond to the two breathers 
(it is easy to check their masses are m t and m s respec- 
tively). The third order pole at 6 = 2ni/3 in S& 2 ,& 2 (#) 
should be interpreted as associated to intermediate vir- 
tual particles in the scattering processes^, as drawn in 
Fig. [21 All other poles in (|3I5[) are redundant and do not 
correspond to additional bound states^. 

The finite-size corrections for a QFT in a large but 
finite volume with periodic boundary conditions are a 
consequence of the vacuum polarization and, when the 
model is integrable, they can be extracted from the ex- 
act scattering data of the infinite-volume theory^. The 
leading corrections to the masses are exponentially small 
and consist of two terms 

Am a (L) = Ami F \L) + Am^(L) 7 (6) 

which roughly can be understood as follows. The first 
term exists in any theory and can be interpreted as vir- 
tual particles "traveling around the world" once before 




FIG. 2: Physical interpretation of third order poles. The 
angle 6 is related to the poles in the intermediate processes 

as 9 = u b,a> + U L' - t 1 "- 

being annihilated again. It has the form 

Amf\L) = - E ^- e " m e b (0)faA6) , (7) 

where f a , b {8) = S a . b (9 + f ) - 1. For large L Eq. Q is 
of order 0(e~ mbL ). The second term in © is present 
only in theories with bound states and is associated to 
virtual processes in which a particle is split into its two 
constituents that "travel around the world" before re- 
combining again into the original particle. If the particle 
a is bound state of two other particles of the theory, say 
6, c, associated to the pole 9 = iitg c , then Ami^L) reads 

b.c 

(8) 

where denotes the Heaviside step function and 

Mabc = m b smu c ab , R abc = -iRese = i u <= ab S a , b (9) . (9) 

As mentioned before the third order pole has to 
be interpreted as formation of virtual particles and 
R abc = Rah'c'Rha'c'Rch'a'- In general, since for small L, 
AE St t(L) ~ 1/L, if the terms above approach the infinite- 
volume limit from below the scaling function will present 
a minimum (the minimum may also appear as a conse- 
quence of the competition of the two terms in 

Since the singlet and the triplet masses are related, 
all the finite-size corrections depend on the scaling vari- 
able m t L (Ref. \l$) . Using formulas il(ilSI) and the 
pole structure of the S-matrix, we find the following ex- 
pression for the leading finite-size corrections, $ s t(£) = 
A^ lt (L)/A^(oo), 

/oo in 
-00 27r 

= V3 (l - 3e- £ / 2 - 36e M / 2 ) , (11) 

with f Q (9) = 2 % /3cosh6»/(2cosh6» - y/3). These scal- 
ing functions are compared to the scaled DMRG data in 
Fig. (one has to replace m t L with the scaling vari- 
able £ = L/£). For each value of 8 the scale factors £ 
(slightly different for triplet and the singlet) and AE t 
have been tuned in order to give the best collapse of 
the curves. Using the numbers found for the triplet, the 
product v = £AEt yields another estimate of the velocity, 
V = 1.19 ± 0.01, in agreement with those in Table [I] 



4 



- f (a) 

- ? 


1 1 1 1 1 


!• 










Singlet - 












Triplet 


I i » i i i i i t ] t 1 t ] 



2 4 6 



10 12 14 



J 2 =-0.03, y,=-0.10 (coinciding) 




Triplet 



A first interesting result of this analysis is that sta- 
ble excitations of the lattice model for J-2 = J2c with a 
dimerization term exhibit a clear scaling behavior over 
the whole range of L for different values of S. Moreover, 
the scaling functions obtained using the QFT description 
capture well the finite-size effects of the spin Hamiltonian 
for I > 4 or 7, depending on which of the two excitations 
we consider. We recall that the analytic results are the 
exact expression of the leading finite-size corrections and 
contain no additional fitting parameters. 

Finally, we have repeated the numerical analysis for 
other values of J-2- Unfortunately in this case the low- 
energy effective field theory is no longer intcgrablc and 
computations of finite-size corrections are much more dif- 
ficult. In principle both the variation of the mass and of 
the S-matrix can be taken into account perturbatively^S. 
The numerical results are reported in Fig. We note 
that the form of the scaling function for the singlet has 
a much stronger dependence on J2 than the triplet one. 
For fixed J2 we find again no appreciable dependence 
on S. 



FIG. 3: (a) Scaling plot of finite-size corrections to the triplet 
and singlet gaps AEt, s (L) at J2 = Ji^. Continuous lines rep- 
resent the analytical predictions 1101 and lilt , while the sym- 
bols refer to DMRG data obtained with the same parameters 
as in Table[I](circles: 5 = 0.01; diamonds: S = 0.02; triangles: 
S = 0.10). The dashed line shows the l/£ behavior at small 
I. (b) Same as in (a) for J2 < J2 C (squares: S = 0.06). The 
singlet data have been rescaled by A_B s (oo) and then shifted 
by +0.4 for clarity reasons. 
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